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We discuss the bosonization of non-relativistic fermions in one space dimension 
in terms of bilocal operators which are naturally related to the generators of W- 
infmity algebra. The resulting system is analogous to the problem of a spin in 
a magnetic field for the group W-infinity. The new dynamical variables turn out 
to be W-infinity group elements valued in the coset W-infinity/if where if is a 
Cartan subalgebra. A classical action with an H gauge invariance is presented. 
This action is three-dimensional. It turns out to be similiar to the action that 
describes the colour degrees of freedom of a Yang- Mills particle in a fixed external 
field. We also discuss the relation of this action with the one we recently arrived 
at in the Euclidean continuation of the theory using different coordinates. 
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Introduction : 

In this paper we discuss the bosonization of non-relativistic fermions in one 
space dimension. This is a very old subject beginning with the classic work of 
Tomonaga [1] in which he succeeded in deriving the equations of motion of the 
collective excitation (sound wave) near the fermi surface, under certain approxima- 
tions. There have been many subsequent studies of this problem [2] . More recently 
there was a revival of interest in this issue in the discussions of the c = l,d = 2 
string theory. As is well known, this model is formulated as a many fermion problem 
in a background potential [3]. The bosonization of this theory has been discussed 
by many authors from different viewpoints [4-11]. The common feature of all these 
formulations is that the basic boson variable is chosen to be local and is essentially 
the bilinear ip^(x, t)ip(x, i). In a recent work [12] we had argued that a more natural 
set of variables for this problem are the generators of the W 7 ^ ^ algebra which are in 
one-to-one correspondence with the "bilocal" operator ip(x,t)ip'(y,t) = 4>(x,y,t). 
This non-local description is also satisfactory from the point of view that this 
model describes a string theory. In the present problem not all of constitute 
an independent complete set of variables. We have shown [12] from an analysis of 
the classical limit that the independent variables belong to the coset G/H where G 
is the ^-group (group of area-preserving diffeomorphisms in two dimensions) 
and if is a Cartan subgroup. As we emphasized in [12], the coset G/H constitutes 
the 'phase space' of the theory. We also presented a phase space action in terms of 
a certain parametrization of the coset variables (see equation (9.18) of ref [12]). In 
the present paper we employ a different parametrization of the coset to discuss the 
bosonization and derive the classical action in terms of these variables. This con- 
stitutes a solution of the bosonization problem. We hope that besides the solution 
of the bosonization problem that we have obtained here this formulation will help 
resolve some of the perplexing questions of the d = 2 string theory, especially those 

* Here and in the following we shall adopt the matrix notation of [12] and not indicate the 
x, y-dcpendcncc explicitly, unless desirable for clarity. Thus <j>(x, y, t) are components of the 
matrix $(t). 
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related to the puzzling differences in the Minkowski and Euclidean formulations 
and of course the issue of black holes. 

It is appropriate at this point to say a few words about the recent developments 
in this subject. The conservation laws of the fermionic model were first presented in 
[6,7]. Their local form and the algebra of the currents was obtained by the present 
authors [13,12] and in [14, 15]. The possible role of the algebra of area-preserving 
diffeomorphisms and higher spin states in two dimensions was pointed out in [16, 
17]. A precise connection of this algebra with the vertex operators of the d = 2 
string theory was discussed in [18,19]. The role of the VF-infinity algebra in the 
matrix model was discussed in [12] and from a different point of view by [14, 21]. 

Gauge theory formulation of the c = 1 matrix model : 

In the gauge theory formulation [12] of the double scaled fermion field theory 
the action is given by 

S[V, ^, A] — J dt V\t)(idt + A(t))*(t) (1) 

where A(t) is a background gauge field and we have used the notation of [12]. For 
the double scaled c = 1 matrix model A xy (t) = 1/2 (d% + x 2 )5(x — y). This action 
has the background gauge invariance 

9(t) - V(t)*(t) 

2) 

A(t) -> V{t)AV\t) + iV{t)d t V\t) 

where V(t) is a unitary operator on the single-particle Hilbert space in which 
ty(t) is viewed as a vector* with components tp(x,t). In ref [12] we identified the 
generators W^ r,s ^ of the symmetry algebra Wi +OQ of the theory. We also pointed 
out that the relevenace of the symmetry algebra for this problem lies in the fact 

* A more appropriate notation might be to write the ket \$?(t) >, which is equivalent to 
thinking of the second-quantized field as a many-body wave-function. In this notation <&(t) 
may be written as \&(t) X ^{t)\. The components <f>(x,y,t) may then be written as 
< af|*(t) >< V(t)\y >= 1>(x,t)tf(y,t). 
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that the generators W^ r,s ^ are themselves the basic (boson) dynamical variables. 
Their definition involved polynomials of the coordinate and momentum of a single 
fermion. Towards the end of ref [12], in an added note, we also introduced a 
generating function for the operators W^ r,s ^> in a "classical" phase space (p, q) by 
the formula 

W(p, q, t) = 1/2 J dx e ipx i) ] {x + q/2, t) ip(x - q/2, t) (3) 

The generators W^ r,s ^s are then given by 

= Jdadb F rs (a, b) W( b -±^, t) (4) 

where 

F rs (a,b) = 2cos(^) (td a r(td b y5^(a,b) (5) 

The functions F rs carry a representation of the W-infinity algebra. 
The operators W(p, q, t) satisfy a simple algebra 

[W(p,q,t),W(p',J,t)]=i sm{^(pq' - qp')) W(p + p' ,q + q' \t) (6) 

In (6) we have explicitly exhibited the string coupling constant g s tr ~ ^, so that in 
the limit /i^Owe regain the structure constants of the algebra of area-preserving 
diffeomorphisms on the plane. The quantum deformation of the classical algebra 
of area-preserving diffeomorphisms has been discussed in [20]. For the present 
problem it was mentioned in [19]. 

Equations (3) and (6) can be understood in a framework that formulates quan- 
tum mechanics in terms of a classical phase space which is due to Wigner [22]. To 
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explain this, let us introduce the Heisenberg-Weyl group with elements 

a _ e -i(qP- P X) 
yp,q — e 



(7) 



where (p, q) is a point in the plane R 2 and P and X satisfy the commutation 
relation [P, X] = —ih. The group multiplication law is given by 

9p«&J = e th(pq '~ qp ' )/2 g P+P > q+q > (8) 
Now noting that g PA = e -*qP/2 e ipX e -iqP/2^ we see that 

W(p, q, t) = 1/2 J dx ^(x, t) g m ^{x, t) (9) 

One can now immediately recognize that the algebra (6) is an immediate conse- 
quence of the group law (8). 

We now point out that the time evolution of W(p, q, t) is governed by a simple 
equation. In ref [12] we had shown that the time evolution of W^ r,s \t) is given by 
the equation 

d t W^ s \t) = e(r - s)W^ s \t) (10) 

with solution W^ r,s \t) = exp(e(r — s) t) W^ r ' s )(0) (e = 1 in Minkowski and —% in 
Euclidean space). An equally simple equation governs W(p,q,t) 

[d t + e(pd q + qd p )]W(p,q,t) = (11) 

We see that the time development of W(p,q,t) is simply a [/(l)-rotation in the 
phase plane. The solution is 

oo 

W(p,q,t)= W m (T) (12) 

m=— oo 

where we have parametrized the phase plane by p = exp(— r) cosh(e#), q = 
exp(— r) sinh(e6>) In ref [12] we had identified r with the Liouville coordinate of 
the continuum string theory. 
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Equations (6) and (11) constitute an exact boson representation of the non- 
relativistic fermion theory under consideration. 

Boson representation in terms of the bilocal operator : 

The discussion in this and the subsequent sections is best done in terms of the 
original bilocal operator $(£). The operator W(p,q,t) is related to it by a fourier 
transform as in (3). The first point to note is that the equation of motion for 
following from (1) implies the following equation for the bilocal operator $(£): 

idtQ(t) + [A(t),$(t)]=0 (13) 

Moreover, the equal time algebra for $(£) following from the anticommutation 
relation of V(t) and ^(t) is 

[$ Q (t), $P(t)] = % J2 Cf^(t) (14) 

7 

where a,(3,--- denote the components (x a , y a ), (xp, yp), ■ ■ ■ and are the struc- 
ture constants 

= -iS(xp - y a )S(x^ - x a )8(y 1 - yp) + iS(x a - yp)8(x 1 - xp)8(y 1 - y a ) (15) 

The sum over 7 in (14) denotes the integration over the variables (x 7 ,y 7 ). The 
structure constants in (15) are related to the standard ones for ^ by a redefi- 
nition of the generators. 

Equations (13) and (14) show that $(t) is like a spin variable and the present 
problem is like that of a spin in a magnetic field. The hamitonian for a quantum 
spin S l in a magnetic field Bj is h = S^Bj and the equation of motion is idfS 1 = 
[h, S] 1 = iC^BjS k using [S^S*- 7 ] = iC^S k . If we define the matrix of operators 
S = S l Ti and a matrix-valued magnetic field B = —B J Tj then the equation of 
motion becomes idfS + [B,S] = in exact analogy with equation (13). In the 
usual case of the quantum spin the TVs correspond to the generators of SU(2); in 
our case the algebra is W^ +(X1 . 
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It immediately follows from (13) that the set of operators tr($(t)) n , n — 1, 2, • • • 
are constants of motion. In fact, these are precisely the Casimir operators for the 
^l+oo algebra (14). This is clear from the fact that tr($(t)) n is invariant under 
$(£) — > V$(t)V~ l where V is any element of the Wjt^-group. Since the above 
rotation is generated by the $ a 's themselves, this implies that 

[tr($(t)) ri ,$ a ] = for all a (16) 

It is the existence of these Casimir invariants that reduces the the number of 
independent variables of the problem. 

The identification of $(t) with a spin-like variable and the constraints coming 
from the Casimir invariants suggest the following parametrization of the bilocal 
operator 

m=U(t)A(t)U-\t) (17) 

where U is an element of the W-jt^ group and A(t) is a diagonal matrix. Since 
tr($(t)) n = tr(A(t)) n , n — 1, 2, • • • are all fixed numbers, it follows that A is a fixed 
time-independent matrix, determined by the values of the Casimir invariants. With 
this ingredient we see that (13) implies the following equation of motion for U{t) 

[K{t),iU- l d t U + U- 1 AU] = Q (18) 

which just states that the constant matrix A commutes with the gauge orbit of 
A. This equation of motion can be derived from the following classical action by 
varying U 

S[U,A,A]= J dttT{(iU~ 1 dtU + U- 1 AU)A} (19) 

A simliar action has previously appeared in the problem of Yang-Mills particle [23] . 
The analogy is exact if the centre of mass motion of the particle is frozen because 
then only the colour degrees of freedom are active. It has also appeared previ- 
ously in the discussions of quantization of baryonic collective excitations around 
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the Skyrme soliton [24,25]. In these works it was shown that for the collective 
excitations the WZW term reduces to the action (19) with A — 0. It is also inter- 
esting to note that there is a natural two-dimensional representation for this action 



It is important to realize that the action (19) is really a 3-dimensional action. 
To emphasize this point we write out the r.h.s. of (19) in detail. Using the notation 
U xy (t) = u(x, y, t),A xy = X(x) 8(x - y), we have 



S[U, A, A] — / dt dx dy X(y) u*(x, y, t) [id t + A(x, d x )\ u(x, y, t) (20) 



An important aspect of the parametrization in (17) is that U{t) is not uniquely 
determined. In fact, U(t) and U(t)g(t) such that g^Ag^ 1 ^) = A are equivalent. 
For generic A, the set of all such elements g{t) forms a Cartan subgroup H of 
^l+oo- The physical degrees of freedom of this system, therefore, belong to the 
coset W^ +OQ /H. One would therefore expect the action in (19) to have the local 
gauge symmetry U(t) — > U{t)g{t). This indeed turns out to be correct [23-25]. We 
also note that (19) has the background gauge invariance of the action (1), i.e (19) 
is invariant under 



U(t) -> V(t)U(t), A(t) -> V(t)A(t)V- l (t) + iV(t)d t V- l (t) (22) 



It is important to realize that the action (19) is first order in time derivative and 
therefore it must be interpreted as a phase space action. Furthermore, as we have 
discussed above, this action has a gauge inviariance which restricts the physical 



[24,25]. 



where u(x,y,t) satisfies the unitarity constraint 




dx'u(x, x', t)u*(y, x', t) = 5(x — y) (21) 



and A = \{d^. + x 2 ) is the background specific to the c = 1 matrix model. 
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degrees of freedom to the coset W^ +00 /H. Thus the action (19) is a phase space 
action with the phase space being the coset. In [12] we arrived at an identical 
conclusion from the analysis of the classical limit, and wrote down a bosonised 
action in terms of a certain parametrization of the coset. As we shall see later, for 
the simple example of SU(2)/U(1) these two actions turn out to be identical. 

How does one determine the constant matrix A appearing in (19)? As we shall 
see below, one can explicitly calculate all the Casimir invariants in the fermion 
field theory described by the action (1). Since the knowledge of all the Casimirs is 
equivalent to the knowledge of A, we can in this way determine A. 

We have proved in [12] that the fermion Fock space in the zero fermion number 
sector forms an irreducible representation of W-jt^ algebra. Therefore, to find out 
the values of the Casimirs it is sufficient to find out their values on the fermi 
vacuum. Consider, for example, the first Casimir C\ = tr $(£). Expanding the 
fermion field in terms of a complete orthonormal basis of energy eigenstates 

${x,t) =^e- iEnt u n {x)c n 

n 

^(x,t)=J2 eiEnt <W C n < 23 ) 
n 

{ c ni c rnl = ^nm 

we can write this Casimir as 

Ci = E c " c n ( 24 ) 

n 

It is clear that C\ on any state measures the total number of unoccupied levels. 
To find its numerical value it is sufficient to find its value on the fermi vacuum. 
Denoting the fermi level by n = no, E no = /i, we see that the value of C\ is 

Ci = £ 1= dE p(E) (25) 

i 

where p(E) is the level density. Similarly, the second Casimir C2 = tr ($(t)) 2 has 
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the expression in terms of the oscillators 



C 2 = Y1 C m4iCn4n (26) 
m,n 

This can be rearranged and put in the form 

C 2 = C 1 (1 + C 1 ) (27) 
where C\ denotes the total number of occupied levels 

C 1 = J2cic n (28) 

n 

On the Fermi vacuum C\ evaluates to 



no ^ 



C 1 = J2 1 = dE P( E ) ( 29 ) 

-°° -oc 

In general, for the n-th Casimir C n = tr ($(t)) n , we find 

C n = C 1 (l + C 1 ) n - 1 (30) 

Equation (30) is ill-defined as it stands, since both C\ and C\ diverge. It is, 
however, sufficient to introduce regulators at the two ends of the spectrum to 
regularize all the Casimirs. The dependence on these regulators should drop out 
of all physical calculations in any well-defined renormalization scheme. 
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Path integral and effective boson field theory : 

So far we have argued in the operator formulation that the action (19) is the 
correct bosonization of the problem. We shall now present arguments within the 
fermionic functional integral for the action (1) to arrive at the same conclusion. 

In the fermion field theory described by the action (1), we are in general inter- 
ested in an arbitrary correlation function of the fermion bilinears tp(x,t)tp^(y,t): 

G(xi,yi,ti] ■ ■ ■ ; x n , y n , t n ) 
= J ^^ t ^(xi,ti)^ t (yi,ti)---^(x n ,t n )^ t (i/ n ,t n )exp(^[^,^ t ,A]) 

(31) 

where Z[A] is the parition function for the action (1). A natural boson representa- 
tion for this correlation function would be the one in terms of the bilocal operator 
(f)(x,y,t). Clearly the classical action for <p(x,y,t) should be such that the corre- 
lation function of (f>(xi,yi, ti), ■ ■ • , <f>(x n , y n , t n ) exactly reproduces (31). One may 
think of this boson action as an effective action for the physical excitations of the 
theory in complete analogy with the mesonic effective actions of QCD. Moreover, 
as we argued in [12], since one can create all the physical states of the theory by 
the action of an arbitrary string of the bilocal operators on the fermi vacuum, such 
a classical bosonic action would represent a complete bosonization of the theory. 

Our strategy for the deriving this boson action is to insert the following identity 
in the fermion field theory partition function which enforces the defintion of the 
bilocal operator <f>(x, y, t) in terms of the fermion fields 

1 = j D$ 6[</>(x, y, t) - 4>(x, t)^(y, t)} (32) 

Using the exponential representation for the delta-functional 

S[(j)(x,y,t)-'ip(x,t)^(y,t)] = J VAexp(i J dt tiA(t)$(t)) exp(i J dt ^ f (t)A(t)^(t))| 

(33) 
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the partition function can be written as 

Z[A] = J V$VAexp(i J dt tr A(t)$(t))Z[A + A] (34) 

where 

Z[A] = J Dtf+PWexp^ J dt^(t)(id t + A(t))V(t) (35) 

Clearly, A is like a source for the fermion bilocal operator in (35), and log Z[A+A] 
is the generating functional of connected correlation functions for these in the 
background A. The integration over A in (34) performs an inverse fourier transform 
of Z[A + A] and hence gives us the desired classical bosonic action. 

Let us now make the following change of variables in the functional integral 
$(t) = U{t)A{t)U{t)-\A{t) -> U(t)A(t)U(t)-\^(t) -> U(t)*(t) (36) 
where A(t) is diagonal. The partition function now becomes 

Z[A] = J VUVA J(A)VAexp(i J dttrA(t)A(t))Z[A + A u ] (37) 

where 

A u = iU^dtU + U~ l AU (38) 

and J (A) is the van der Monde Jacobian coming from the change of variables from 
the hermitian matrix $ to its eigenvalues A and "angles" U. We do not expect any 
anomalies in the fermionic measure because we are dealing with non-relativistic 
fermions. 

To be more precise the integration over U should really be over a coset as is 
clear from the counting of the degrees of freedom of $ on the one hand (remember 
it is a hermitian matrix) and A and U on the other hand. But, as we have men- 
tioned above, the action of the U variables turns out to be (19) which has a gauge 
symmetry which enables us to gauge away the extra degrees of freedom. Therefore 
the integration over the C/'s can be extended over the whole group. 
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We emphasize that the measure for [/-integration is the usual Haar measure 
and that there is no [/-dependent contribution from the Jacobian. As a result the 
entire classical action for U comes from the coupling of the dummy gauge field 
A(t) to the eigenvalues A. In fact, this action, obtained by simply shifting in 
(37) the integration variable A by A u , is identical to (19) provided we can show 
that A is independent of time. We have seen that this is true in the operator 
formalism where we argued that A is determined in terms of the Casimirs which 
are time-independent numbers. In the functional integral, integration over the 
dummy gauge field A and the fermions gives rise to the delta functional 5[c^A]. 
Let us explain in some detail how that happens. 

Since A is diagonal it couples only to the diagonal elements of A. So the only 
source of dependence on A in the A -\I> integration in (37) is that over these diagonal 
elements A xx (t) of A. This problem is like QED with the gauge field coupled to an 
external source J M and without the kinetic term for the gauge field. If the external 
source has a longitudinal piece then integration over the gauge degrees of freedom 
of the gauge field does not give rise to the gauge volume but to a delta functional 
of d^J^. In the present case, we can decompose A xx (t) into gauge and physical 
degrees of freedom as 

A xx (t) = d t e{x,t) +a(x) (39) 
The A-ty functional integral in (37) can then be written as 

J Ve(x,t) J Va{x) exp(i J dt dx (a(x)A(x,t) - e(x,t)d t A(x,t)) Z[a,e] (40) 

where 

Z[a,e] = JvaJ vtfVVexp(i J dt \t){id t + d t e{t) + a + A(t))*(t)) (41) 

In (41) e(t) and a are diagonal matrices with components e(x,t) and a(x), and the 
gauge field A(t) has only off-diagonal elements non-zero which are being integrated 
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over. Since e(t) is diagonal we can remove it from the r.h.s. of (41) by the change 
of variables 



which does not affect a since it is also diagonal. In effect this means that Z[a, e] = 
Z[a] is actually independent of e. The e integration in (40) is, therefore, trivial and 
gives rise to the delta-functional S[dtA(x,t)] as stated above. 

We now come to the integration over the constant matrices A in (37). That 
there is an integration over A may seem surprising in view of our demonstration 
earlier that A is a fixed c-number matrix determined by the values of Casimirs. 
What this means is that in the functional integral there must be delta-functions 
which enforce this. We must remember that the classical action (19) is a phase 
space action and therefore physical states are functionals of only "half the coset 
variables (the coordinates). Our demonstration that the Casimirs have fixed values 
on such wave-functionals implies that these delta-function constraints arise in the 
functional integral only after integrating out the other "half of the variables (the 
momenta). 

We may therefore write the partition function in (37) as 



and it is understood that the A-integration will eventually get frozen. We can 
similarly argue that the boson representation for the correlation function in (31) is 



A(t) -> e ie ^A(t)e- l <^ 



(42) 




(43) 



where 




(44) 



G(x\,yi,ti; ■ ■ • ;x n ,y n ,t n ) 




-iS[U,A,A] 



(45) 
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Since the A-integration in both the numerator and the denominator in (45) is 
eventually frozen, we may effectively write the above as 



G{xi, yi, h; •••;*„, y n , t n ) = ______ 1 

(46) 

This represents a complete bosonization of the theory and we have thus derived 
the effective classical bosonic action (19) in the functional formulation. 
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The example of SU(2)/U(1) 

Finally, an important question to ask is how to paramatrize U in terms of 
physical coordinates and momenta. This requires fixing a gauge to restrict U to 
the appropriate coset and choosing a polarization. A simple illustrative example 
is that of SU{2) /U(l) [23, 26]. For the simple choice of A = A = 03, the action 
(19) in this case becomes 

S = J dt tr {iU^dtUaz + U'^Uaz) (47) 

where U (t) are SU(2) group elements. It is easy to check that this action has the 
local gauge symmetry U(t) — > U(t) exp(i9(t)as). Therefore the phase space in this 
case is the manifold SU{2)/U{\) = S 2 . There is a natural action of SU(2) on S 2 
if we think of S 2 as CP 1 . In other words let us represent points on the sphere 
by a pair of complex numbers (21,22) with the understanding that (21,22) and 
(A2i,A22) represent the same point (A is any non-zero complex number). Except 
when 22 = we can describe the point (21, 22) by the "inhomogeneous coordinate" 
2 = 21/22. Now an SU{2) matrix acts on the pair (21, 22) by sending it to (021 + 
bz2, C21, d,Z2), ad — be = 1. This implies that 



ax ~\~ b 

U:z^ -, ad-bc=l (48) 

cz + d v ' 



We would now like to parametrize points of the sphere by group elements in the 
following fashion. Let us fix a point on the sphere, say (21,22) = (0, 1) or equiv- 
alently 2 = (we can think of this point as the North pole if we identify this 2 
with the stereographic coordinate). Let us ask the question which group element 
U takes us from the North pole to a given point 2. Of course the answer is not 
unique because any digaonal matrix (6 = c = 0) will keep 2 = at 2 = 0, according 
to (48). This is simply a reflection of the fact that the 2-sphere is a coset of SU{2). 
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Let us choose the following representative element 




U(z) 



Clearly this sends z = to z. 

Once we make this gauge choice it is simple to work out the action (47) in 



terms of z,z. 



Uzd t z - zd t z) + 1 - 2zz , 

S = 2 dt ^— -t— 49 

1 1 + zz K J 



In terms of the variable w = z/y/1 + zz the action becomes 

S = 2 J dt [^(wdtw - wd t w) + 1 - 3ww] (50) 

This may be compared with the action we had written down in [12] (see eqation 
(9.18) and the footnote following it). These two actions agree, apart from a trivial 
shift of the hamiltonian by a constant number and a rescaling of time. 

Concluding Remarks 

In summary we have presented a classical action for the bosonic excitations of 
the c = 1 fermion field theory. The action is a phase space action written in terms 
of VF^oq group element U(t) valued in a coset W^ +00 /H. The action is naturally 
written in three dimensions. An important open problem that remains is to find 
a parametrization of the coset that lends itself to a two-dimensional spacetime 
interpretation of the spectrum of the model. To this end, it might be appropriate 
to use "mixed" bases such that u n (x,t) =< x\U(t)\n >, u^(x,t) =< n\U*(t)\x >, 
where \n > are the eigenstates of A, A\n >= E n \n>. Taking A to be diagonal in 
the \n > basis, i.e. A = J2 n \ n \n >< n\, the action reads in the mixed bases 

S = ^2^ n J dt dx u^(x,t)[id t + A(x,d x )]u n (x,t) (51) 
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The unitarity constraints on the u n (x,t) are 

u n( x , t)un(y, t) =5(x - y) 

, U ^ (52) 
/ dx u n [x, t)u m [x, t) =S nm 

Using (52) one can check that for generic A n (51) has the gauge invariance u n (x, i) — >| 
exp(i0 n (t))u n (x, t). For non-generic \ n (i.e. if two or more eigenvalues coincide) 
the gauge-invariance is enhanced. It is remarkable that the equation of motion 
for the bosonic variable u n (x,t) is the same as that satisfied by the fermion field. 
However, here the boson variables are constrained by (52). 

A detailed analysis of the model along the above lines is in progress and will 
appear elsewhere. 

Acknowledgement : One of us (S.R.W) would like to thank I. Bakas for discus- 
sions. 
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